Journal of the Mechanical Behavior of Biomedical Materials

The human masseter muscle is one of the primary muscles responsible for mastication and mandibular movement; however, its intrinsic mechanical properties remain insufficiently characterized. In this experimental study, the nonlinear, viscoelastic, and history-dependent behaviour of the human masseter muscle was investigated using ex vivo uniaxial cyclic tensile tests. The masseter muscle samples prepared from fresh and formalin-preserved cadavers were tested under two loading protocols: a continuous stretch protocol with increasing stretch levels and a constant stretch protocol with repeated loading to a fixed maximum stretch. Tests were conducted at two strain rates, and their influences on the mechanical behaviour of the tissue were examined. The effect of formalin preservation was also investigated. The results showed that the stiffness of the tissue increases for formalin-preserved samples. Under cyclic loading, the features including energy dissipation, stress-softening, residual deformation, and cyclic conditioning progressively changed during the initial loading cycles and reached stabilization during the final cycle. These findings provide experimental evidence that the human masseter muscle exhibits nonlinear, viscoelastic, and history-dependent mechanical behaviour under cyclic tensile loading. The experimental data obtained in this study may be used for biomechanical modelling of the human masticatory system and the development of constitutive models for cranio-maxillofacial surgical simulation, prosthetic design, and facial soft-tissue biomechanics. Statement of significanceThe masseter muscle is one of the primary muscles of mastication. To address the current gap in craniofacial biomechanics that has largely focused on the mechanical characterization of the masseter muscle based on imaging techniques or monotonic loading, this study quantifies the nonlinear and viscoelastic mechanical response of masseter tissue under cyclic continuous and constant stretch loading, including strain-rate and preservation effects. The results show that the mechanical behaviour of the masseter muscle, including stiffness, hysteresis, stress-softening, and residual strain behaviour, is strongly influenced by strain-rate and formalin preservation. The experimental results provide mechanical data for constitutive modelling of the masticatory system with applications in cranio-maxillofacial surgical simulation, prosthetic design, and facial soft tissue modelling.

The macroscopic biomechanical characteristics of soft and ultrasoft tissues, such as articular cartilage and vocal folds, significantly determine their physiological function. Treatments of widespread tissue degradations due to osteoarthritis in the knee or vocal fold impairment remain an unresolved challenge. For the design of implants for tissue repair after injury or disease, it is key to thoroughly understand the unique biomechanical properties of native tissues and potential substitute materials. We use multimodal mechanical testing methods combined with hyperelastic nonlinear continuum mechanics modeling, and finite element simulations to determine the macroscopic behavior of surrogate materials for human articular cartilage in the knee and human vocal folds. Our cyclic loading experiments reveal qualitative similarities for both tissues and their surrogates, including a nonlinear stress-strain behavior, hysteresis, and conditioning. We demonstrate the tunability of biomimetic and biosimilar stiffnesses of synthetic articular cartilage and vocal fold surrogates through tissue-specific process-material combinations. Our results demonstrate the feasibility of synthetic metamaterials in replicating essential passive biomechanical functions with great potential for future treatment options.

Acute subdural hematoma (ASDH) is a prevalent injury with high mortality and morbidity, often resulting from bridging vein (BV) disruption secondary to cortical relative motion. As a thin membrane enveloping the brain surface and anchoring BVs, the pia mater is hypothesized to play a critical mechanical role in cortical response and hence ASDH pathogenesis. Finite element (FE) head models are valuable tools to predict ASDH occurrence during impacts. However, the pia mater is often represented as an elastic material in existing FE head models, despite experimental evidence reporting its nonlinear mechanical behavior. In this study, both linear (Young's modulus of 11.5 MPa) and nonlinear (the stress-strain curve derived from pial tension tests) material models of the pia mater were implemented in one FE head model. The models were subjected to three experimental impact loadings, one of which was known to cause ASDH and two of which were not. Results demonstrated that, across all simulated impacts, the model with nonlinear pia mater properties predicted larger cortical displacements and BV responses than the linear model. For the impact with known ASDH occurrence, the predicted BV strain was 0.17 for the nonlinear model and 0.094 for the linear model, with only the former approaching the reported rupture strain range of the BV-superior sagittal sinus complex (0.29 {+/-} 0.13). These findings verified the mechanical importance of the pia mater in cortical responses and hence the prediction of ASDH, suggesting that conventional linear pia modeling might over-constrain cortical motion, leading to underestimation of BV strain and ASDH risk. The current study supported the adoption of experimentally derived nonlinear pia mater properties in FE head models to improve the reliability of ASDH prediction.

Purpose: Fibrotic remodeling of the uterus, associated with aging, disease, and environmental exposures, alters collagen organization and tissue stiffness, yet how these changes influence organ-level mechanical behavior remains poorly understood. Glutaraldehyde (GA)-induced collagen crosslinking was used as a controlled surrogate for fibrotic remodeling to determine whether image-informed inverse finite element analysis (iFEA), combined with inflation testing and micro-computed tomography (microCT), could detect and quantify the resulting changes in uterine constitutive behavior. Methods: Murine uteri (n = 6 untreated, n = 6 GA-crosslinked) underwent volume-controlled balloon inflation with simultaneous microCT imaging to quantify deformation of the inner and outer wall boundaries for iFEA. Specimen-specific Gasser-Ogden-Holzapfel (GOH) finite element models were optimized by adjusting model parameters to reproduce experimentally measured wall contours throughout inflation. Model performance was evaluated using contour root mean square error (RMSE), and parameter identifiability was assessed through sensitivity analyses. Results: GA treatment significantly increased inflation work, linear stiffness, and maximum inflation resistance (p < 0.001). The iFEA framework accurately reproduced experimental deformation (RMSE < 3%) and revealed significant increases in the estimated GOH parameters C10 (9.2-fold), k1 (2.0-fold), and k2 (2.7-fold), consistent with increased effective tissue stiffness and a shift toward earlier collagen fiber recruitment. Sensitivity analyses demonstrated unique, well-defined minima for all parameter combinations. Conclusion: Image-informed iFEA provides a quantitative framework for relating collagen remodeling to organ-level uterine mechanics through specimen-specific constitutive parameter estimation. This approach establishes a foundation for investigating the mechanical consequences of uterine fibrosis and other remodeling processes.

Spinal cord injury (SCI) causes rapid and severe bone loss in the paralysed lower limbs, particularly at the distal femur and proximal tibia, where fragility fracture risk is high. In vitro nanoscale vibration at 1 kHz has been shown to promote osteogenic differentiation and inhibit osteoclastogenesis, suggesting potential as a targeted mechanical intervention. This study aimed to develop and evaluate a wearable device for delivering and monitoring localised nanovibration at the distal femur in individuals with SCI. The device delivered continuous sinusoidal nanoscale stimulation at 1 kHz via a bone-conduction transducer, with an opposing accelerometer used to monitor transmitted vibration in real time. Design and target-site selection were refined through two healthy-volunteer investigations comparing the distal femur, proximal tibia, and distal tibia. Bovine femur experiments characterised vibration transmission under controlled benchtop conditions. Preliminary repeated-use feasibility was assessed in one individual with motor-complete SCI. Healthy volunteer testing showed that although the ankle initially produced the highest transmitted amplitudes, these were highly variable, and positioning was inconsistent. Within the knee region, the distal femur provided the most practical and repeatable site for a wearable application. In bovine femur experiments, scanning laser vibrometry demonstrated measurable vibration on the condylar surface opposite the transducer, and depth-resolved measurements confirmed that nanoscale vibration remained detectable within bone. A gel interface layer reduced the transmitted amplitude. In the feasibility evaluation, 61 sessions were completed over 14 weeks, with logged accelerometry confirming repeated nanoscale vibration transmission. These findings establish feasibility and support further device optimisation and translational studies.

The cytoskeleton is a dynamic biopolymer network whose shear rheological properties are crucial for cellular physiology and pathology. However, its mechanical behavior spans multiple spatiotemporal scales, and the coupling of dynamic remodeling and viscoelastic dissipation mechanisms poses a challenge for traditional models to comprehensively capture complex cellular responses. This study aims to establish a multiscale cytoskeletal network model that integrates the bio-chemo-mechanical properties of local linked proteins, the viscoelasticity of actin filaments, and their deformation states. Developing a boundary-modified finite element method with an incremental iterative algorithm, we demonstrated the dynamic remodeling of network and the resultant rheological properties of cytoskeleton by extending the predictive time scale to one thousand seconds. The results not only reproduced the short- and intermediate-term power-law creep behavior and long-term strain plateau response of the cytoskeletal network observed in shear rheological experiments, but also indicate that the synergy among the chemo-mechanical coupling of cross-linked proteins and the bending-to-tension transition of actin filaments govern both the network remodeling and its power-law response evolutionary, whereas the steady-state properties of actin filaments determine the long-term network behavior. Simulations of cancerous and drug effects show that cancer-induced softening and reduced filament viscosity lead to accelerated cytoskeletal responses and decreased apparent shear modulus, respectively; and drug-enhanced filament prestress, along with promoting association or inhibiting dissociation of cross-linked proteins, can effectively increase the steady-state shear modulus. These findings advance the understanding of the spatiotemporal evolution and pathological mechanisms of cellular mechanical responses and provide insights for regulating polymer network performance.

PurposeComputational models of the uterus during pregnancy enable analysis of electro-chemo-mechanical pathways to predict labor timing and guide treatment planning. We aim to develop a robust image-based modeling pipeline to investigate uterine passive mechanics during late pregnancy. MethodsA parametric model of the uterus and cervix was created using a patients MRI measurements at 38 weeks of gestation. Inspired by advances in cardiac mechanics models, we created Laplace-Dirichlet solutions to inform tissue domains, fiber structure within the uterus and cervix, and spatially varying Robin boundary conditions. Prior imaging and mechanical testing data were used to fit material parameters. Boundary condition parameters were tuned to match the displacements of a previously established approach that employed contact with surrounding tissue. The tissue mechanical response to a physiologic load was assessed across varying material properties and fiber architectures. ResultsDiscrepancies in nodal displacements between the current approach and the contact-based model were limited to 3.4 {+/-} 1.8 mm, yielding nearly 90 % computational savings. Uterine tensile strains were more sensitive to ground substance elastic modulus (E) compared to fiber properties. Reduced E and fiber stiffness increased cervical strains and compression. Fiber dispersion and architecture modulated the opening of the cervical internal ostium but had a reduced impact on compression. ConclusionWe developed a novel workflow for modeling passive uterine mechanics, informed by patient-specific measurements and in vitro mechanical tests. The robust workflow may prove useful for studying labor progression and conducting longitudinal studies to enhance our understanding of normal and pathological pregnancies.

Successful instrumented fusion of the lumbar spine is a complex surgical challenge, with positive patient outcomes dependent on careful surgical planning. Material selection is of critical importance to a mechanical construct supporting successful spinal fusion. Therefore, the aims of this study were to (a) evaluate the potential clinical use of finite element analysis (FEA) and (b) conduct a retrospective mechanical analysis of different implant materials in patients having undergone spinal fusion using FEA. Our methodology involved segmenting the spine from post-operative computed tomography (CT) image data from patients with previous spinal fusion. FEA models representing post-surgery cases were developed and different biomechanical loading conditions such as compression, flexion, bending and extension whilst testing pedicle screws of different materials were simulated. Patient specific finite element models were created, and biomechanical analysis were completed for all three patients. Polyetheretherketone (PEEK) constructs typically demonstrated lower peak implant stress when compared to titanium constructs for all spinal fusion levels. Furthermore, increasing the spinal fusion level resulted in significant differences in the maximum von Mises stress within both the bone and the instrumentation, whereas the 2-level fusion exhibited comparable stress levels in the bone irrespective of the instrumentation material. This pilot explores the potential of FEA as a clinical tool for assessing device and bone stresses. In our cohort, different materials can influence the stresses in both the instrumentation and the instrumented vertebrae, suggesting FEA can be useful pre- operative tool with regards to instrument selection and post-operatively to assess instrumentation and bone stresses. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=74 SRC="FIGDIR/small/742715v1_ufig1.gif" ALT="Figure 1"> View larger version (33K): org.highwire.dtl.DTLVardef@55959aorg.highwire.dtl.DTLVardef@d0b9d6org.highwire.dtl.DTLVardef@158c348org.highwire.dtl.DTLVardef@7ce828_HPS_FORMAT_FIGEXP M_FIG C_FIG

Repetitive low-intensity blast (LIB) exposure has been identified as a probable cause of mild blast-induced traumatic brain injury (mbTBI) and a chronic injury risk to U.S. military personnel. However, the human brains biomechanical response to this loading regime remains poorly characterized. Using blast-exposure-validated 3D models of human anatomy, we simulated the intracranial tissue response to blast pressure typically experienced by Warfighters during weapons training. Two scenarios were evaluated, a single-dose exposure and a repetitive-dose exposure, to study intracranial pressure (ICP), shear strains, and spectral content. Ansys LS-DYNA was used to generate planar blast waves with peak overpressures of 4-90 kPa and positive phase durations of 2.2-10 ms. Single exposures produced ICP ranging from 4.7-112.7 kPa, dependent on dose and positive phase duration. Under repetitive LIB exposure, peak ICP increased by 8-26% relative to single exposures, with an increase of high-frequency components (>2 kHz). These results demonstrate that LIB can produce measurable intracranial responses that are amplified through repetition, producing pronounced spectral content and elevated pressures despite low strain levels. This study underscores the need to further investigate cumulative dose effects and the value of computational approaches to clarify hypothesized mbTBI mechanisms in operationally relevant conditions.

Purpose: Quantifying aortic morphology is central to surgical planning for thoracic endovascular aortic repair (TEVAR), yet no consensus exists on how best to represent three-dimensional aortic shape for outcome prediction. Two broad strategies have emerged: statistical shape analysis (SSA), which relies on statistical methods and dimensionality reduction to capture the most significant shape modes, and geometrically-informed approaches that extract descriptors grounded in differential geometry. Here, we directly compare these paradigms on a cohort of 290 CTA scans classified by surgical outcome (non-pathological, successful TEVAR, failed TEVAR). Methods: For the geometrically-informed approach, we use a two-dimensional feature space using normalized fluctuation in integrated Gaussian curvature $\widetilde{\delta K}$ and mean aortic radius $R$. For SSA, we construct a point-cloud shape model with dimensionality reduction using Principal Component Analysis (PCA) and evaluate classification performance as a function of the number of retained principal components. Results: SSA's leading principal components encode variations in global aortic size and are statistically redundant with ($R$, $\widetilde{\delta K}$), yet they lack a one-to-one correspondence with interpretable anatomical quantities. Testing on an unseen, independent dataset reveals that the geometrically-informed approach provided better generalizability than SSA. Using Gaussian process classification with 10-fold cross-validation, we find that the geometrically-informed approach achieves a higher weighted $F_1$ score than SSA achieves with up to 20 principal components. While SSA's full-dataset accuracy rises above 90\% with increasing dimensionality, this gain is driven by overfitting rather than genuine discriminative power. Conclusion: These results demonstrate that geometrically-informed descriptors offer a more interpretable, robust, and clinically translatable framework for aortic disease classification than data-driven statistical shape representations.

Silicones are elastomers that have a wide variety of uses, including biomedical applications such as the coating of biomedical devices and as implants. Soft silicones with mechanical properties similar to those of biological tissues have particularly gained use as substrates for cell culture in mechanobiology studies. In this context, it would be desirable to be able to alter their surface mechanical properties with a relatively simple physical treatment. While deep ultraviolet (deep UV) or ultraviolet C (UV-C) treatment has been previously used as a surface treatment method for stiffer silicones formulations, the effect of this treatment on soft silicones relevant for mechanobiology applications is still uncharacterized. We first used nanoindentation to determine the Youngs modulus of two types of soft silicones, Qgel and GEL-8100/Syl (GEL-8100 with Sylgard-184 crosslinker), both with initial moduli in the kilopascal range. We show that nanoindentation in the presence of 1% sodium dodecyl sulfate avoids adhesion between the nanoindentation glass probe and the soft silicones. After deep UV exposure in the presence of air, nanoindentation revealed that the apparent Youngs moduli of the soft silicones Qgel and GEL-8100/Syl increased by 70% and 33%, respectively. The bulk rheology of the soft silicones were not affected, suggesting that this corresponds to a surface stiffening effect with a topical stiffening of at least several hundred kilopascals. Energy-dispersive X-ray spectroscopy results show an increase in the mole fraction of oxygen, consistent with oxidation of the surface. Attenuated Total Reflectance Fourier-Transform Infrared spectra show evidence of Si-OH group formation in GEL-8100/Syl and silicon sub-oxide formation in Qgel. Consistent with this, water contact angle measurements show enhanced hydrophilicity after deep UV treatment. Our results have implications for using soft silicones as substrates in mechanobiology studies and in processes where deep UV light is used in the surface treatment of soft silicones.

Introduction: Metastatic spinal lesions substantially alter vertebral mechanical properties and increase fracture risk. Computed tomography (CT) based finite element (FE) models can estimate vertebral strength, but their accuracy depends on how CT derived material properties are represented. This study evaluated the effect of two material grouping strategies on simulated strength and stiffness in metastatic vertebrae. Methods: We compared Adaptive Clustering (AC) with Uniform fixed width grouping in 44 vertebrae from 11 donors (8 osteolytic, 12 osteoblastic, 12 mixed, 12 no observed lesion (NOL)). FE models were generated based on CT scans with 2 to 500 material groups and compared for material mapping error and simulated strength and stiffness. Overall and lesion stratified agreement with experimental measurements was assessed in an exploratory analysis. Results: AC showed significantly lower Young's modulus root mean square error than Uniform (p < 0.05). Simulated strength and stiffness stabilised by 50 material groups. At 50 groups, simulated strength showed moderate correlation with experimental strength overall (R2 = 0.57), strongest in NOL vertebrae (R2 = 0.82) and lower in lesion-bearing vertebrae (R2 = 0.4-0.59). Stiffness showed weaker correlation overall (R2 = 0.27), highest in NOL vertebrae (R2 = 0.48) and negligible in mixed lesions (R2 = 0.007). Bland Altman analyses indicated systematic underestimation of experimental fracture load. Discussion: AC improved material-mapping fidelity, whereas increasing material groups beyond 50 had little influence on simulated strength or stiffness. Numerical stabilisation therefore did not imply experimental accuracy. Lesion stratified findings were exploratory and should be interpreted cautiously because of limited subgroup sizes.

In the current clinical practice, the diagnosis of spinal disorders and their surgical planning are critically based on imaging data. To complement this data, patient-specific finite element models have been developed and showed to be powerful tools for evaluating spine mechanics. Most of them rely on Computational Tomography (CT) scans - which have a high resolution but are seldom available in routine clinical practice - while only a recent few models are on less invasive Magnetic Resonance Imaging (MRI). Yet, despite the proliferation of these computational models, encompassing detailed anatomical and functional information, the rheological assumptions they are built upon are based on tissue-sample mechanical response data, which leaves a gap in the quantitative analysis on how such assumptions influence the macroscopic response of a functional spinal unit. Aiming at addressing these shortcomings, the main purpose of this work is to introduce a quantitative computational assessment of the macroscopic impact of commonly adopted rheological models - from linear elasticity to fiber-reinforced nonlinear hyperelasticity - in several loading conditions, focusing on a lumbar unit which is considered as a typical benchmark system. We also propose a reconstruction procedure to accurately describe subject-specific anatomy from MRI data, including the intervertebral disc and its nucleus pulposus. Bones are modeled as linear elastic media, whereas for the AF, we consider three different mechanical models - namely, isotropic linear elasticity and the Holzapfel-Gasser-Ogden model with and without fiber reinforcement. Model verification on an idealized geometry demonstrates numerical consistency, while parametric orthostatic simulations highlight the need for nonlinear formulations to capture anisotropy and strain-stiffening behavior of the intervertebral disc. Then, we carry out flexion, lateral bending, and torsion tests on a subject-specific reconstructed functional unit, for which we provide parametric analysis in terms of momentum magnitude and resulting range of motion. These tests further confirm the need for a nonlinear rheology of the annulus fibrosus and provide a quantitative assessment of the differences between the constitutive laws considered. Moreover, successful comparisons with the literature, in terms of macroscopic deformation under several loading conditions, serve as partial validation for our computational model.

Mechanical forces have a major effect on cell behavior. Most cells in vitro are grown under static conditions on hard tissue culture plastic, conditions that do not accurately reflect living tissues. The ability of cells to sense and respond to mechanical forces is essential for key biological processes, including development, proliferation, and migration. Disruption of the ability to respond to mechanical forces are known to be a critical factor in many diseases, including cardiovascular disease, progeria, and cancer. Here, we present the design, fabrication, and biological testing of a custom-built cell-stretching device that applies controlled biaxial strain to cells cultured on a polydimethylsiloxane (PDMS) membrane. We then used this device to examine how cells respond to strain. In response to biaxial strain, MCF-7 cells activated the mechanosensitive immediate early gene (IEX-1), with its expression increasing significantly after 1 and 3 hours of stretching. Cells exposed to mechanical strain also remodeled their cytoskeleton in a direction-dependent manner. Under uniaxial strain, actin filaments reoriented perpendicular to the stretch direction, whereas biaxially stretched cells do not promote directional reorientation, but instead appear to reinforce actin at the cell periphery. Similarly, cells under uniaxial strain exhibited changes in nuclear orientation and shape that were not observed under biaxial strain. Nuclear area remained unchanged in either strain condition. These results highlight that the biaxial stretcher can be used to apply strain to cells, and that cells respond differently to biaxial strain compared to what has been reported for uniaxial strain.

Background and objectiveMusculoskeletal joint disorders often show inconsistent relationships between structural degeneration and nociceptive pain. Temporomandibular joint (TMJ) disc displacement represents a clinically relevant model for investigating the structure-function-pain relationship. This study aimed to develop a multiscale computational framework integrating biomechanics, three-dimensional (3D) neural morphology, and electrophysiology to quantitatively link TMJ structural alterations, biomechanical loading, and peripheral nociceptive activation. MethodsStrain distributions in the TMJ disc and retrodiscal tissue during mouth opening and clenching were computed in ArtiSynth under varying degrees of displacement. Human TMJ 3D nerve architecture was reconstructed using porcine TMJ nerve mapping data as an anatomical reference, and ion channel dynamics were implemented in NEURON. Model coupling was achieved by mapping biomechanical strain fields onto nociceptor membranes to simulate mechanosensitive currents and action potential propagation to the trigeminal ganglion. ResultsAnterior DDwoR induced a severity dependent strain pattern in the TMJ disc and retrodiscal tissue, including posterior redistribution, increased strain magnitude, prolonged activation, and broader retrodiscal tissue involvement. Displacements of 4, 6, and 8 mm produced larger mechanosensitive currents, broader terminal depolarization, and higher trigeminal firing rates during mouth opening (6, 18, and 28 Hz) and clenching (8, 20, and 28 Hz), whereas 0- and 2-mm displacements produced negligible neural activation. ConclusionsThis study establishes a multiscale biomechanical-electrophysiological framework linking TMJ structural alterations to peripheral nociceptive activation. The framework quantitatively connects macroscale strain patterns with microscale neural activation, suggesting that anterior disc displacement may amplify peripheral nociceptive signaling by increasing the overlap between elevated strain and densely innervated retrodiscal tissue.

Snake venom composition and its contribution to toxic effects has been heavily researched, but there is a comparative lack of information on venoms fluidic properties and their relationship with fang morphology during the envenomation process. Understanding how venom flows through a fang can shed light on bite site dynamics and potentially explain bite symptoms. In this article we first conduct a broad comparative test of the rheological properties of venom from thirteen snake species, including multiple viperid and elapid snake species, revealing a shear-thinning non-Newtonian flow behavior in all studied species. However, we have not observed strong phylogenetic signal in venom fluidic properties, suggesting that flow properties may vary independently of evolutionary relationships between snake species. Second, we demonstrate that snake venoms fluidic properties can be modeled by other inexpensive, safe, and abundant shear-thinning surrogate fluids. We found that aqueous solutions of bovine serum albumin protein and xanthan gum are useful venom mimics, matching the rheological behavior of venoms from the studied snake species across a range of relevant shear rates. We further evaluated the performance of these snake venom mimics in a simulated venom delivery system, showing good and robust mimetic control of the flow properties as a function of applied pressure. By elucidating the fluidic properties of snake venom and providing a non-toxic, scalable surrogate fluid model to be used in further studies, we provide the biomedical, toxicology, evolutionary biology communities with a tool to study envenomation physics in an inexpensive and safe fashion. We suggest it is possible to design species-specific venom mimics that facilitate research on the biomechanics and fluid dynamics of venom delivery via snake bites, and inform the design of bioinspired puncture and injection devices.

Keratin proteins are fundamental structural components of hair fibers, contributing to their mechanical resilience, elasticity, and fracture resistance. However, systematic molecular-scale characterization of keratin unfolding mechanics across protein types remains limited, restricting the connection between protein-level deformation mechanisms and hierarchical hair fiber mechanics. Here, we establish a comparative molecular-dynamics-based framework for characterizing the unfolding behavior and nanomechanical response of a curated dataset of 51 keratin proteins. We conduct implicit atomistic molecular dynamics (MD) simulations, including equilibration and steered molecular dynamics (SMD) under four accelerated pulling velocities, to quantify unfolding forces, energy absorption, and structure-property relationships. These accelerated pulling conditions are interpreted as computational probes of relative molecular-scale trends, rather than direct reproductions of experimental hair-fiber strain-rate regimes. Across these accelerated SMD conditions, the simulations show rate-sensitive increases in unfolding force and energy absorption, consistent with constrained molecular relaxation during faster molecular pulling. Stronger correlations between nanomechanical properties and molecular descriptors emerge at higher pulling rates, and the nanomechanical responses of different keratin types (Type I and II) are also compared. The findings provide molecular-level insights into protein unfolding mechanisms that may contribute to the mechanical behavior of hierarchical keratin structures. This study establishes a quantitative framework for comparative keratin unfolding mechanics, providing molecular-level descriptors for future multiscale modeling of hair fiber behavior. These results support applications in biomaterial design, hair fiber durability analysis, and bioinspired material engineering. Future work will integrate these nanomechanical descriptors with fiber-level mechanics and machine learning-based keratin design.

Identification and treatment of pregnancies at risk for preterm birth is a central challenge in obstetric research. Many of the known causes of preterm birth originate from mechanical failure in reproductive tissues. To better understand the biomechanical environment of the gravid uterus and its potential contribution to preterm birth, this computational study presents a parametric method for modeling maternal reproductive anatomy during the early second trimester. A finite element modeling approach was built using existing sonographic measurements from early second-trimester maternal anatomy and material properties from published mechanical tests. We applied the same physiologically relevant intrauterine pressure to all models and quantified the resulting tissue stretch. The sensitivity of the stretch in the proximal cervix was explored by varying material properties and sonographic maternal anatomy dimensions. Cervical material properties, particularly the fiber stiffness modulus and ground substance Youngs modulus, were found to have the greatest effect on proximal cervix stretch compared to other material properties and sonographic dimensions. Among the sonographic dimension measurements, those defining the region surrounding the proximal cervix had the greatest effect on proximal cervix stretch, including the curvature of the posterior uterine wall and the thickness of the lower uterine segment. The computational modeling approach presented here enables future patient-specific studies of gravid reproductive tissues to elucidate differences between individuals who do and do not deliver preterm. Additionally, this study is foundational for building digital twins to support future virtual clinical studies on diagnostic and therapeutic device design to prevent preterm birth.

The incidence of vascular complications differs among dermal filler formulations. This study aimed to evaluate the embolic risk associated with Neofilera(R), a filler composed of carboxymethyl cellulose and polylactic acid microspheres, following intra-arterial injection. The central auricular arteries of rabbits were injected with Neofilera(R) at volumes of 0.1 mL or 0.2 mL under various conditions: normal saline (control) and Neofilera(R) diluted at ratios of 1:5 (Group 1), 1:10 (Group 2), and 1:15 (Group 3). The presence of transparent emboli was assessed immediately after injection, while tissue necrosis (percentage and area) and histopathological alterations were evaluated on days 1 and 7 post-injection. Relative to controls, Neofilera(R) administered at 0.1 mL dispersed within minutes and did not induce significant tissue necrosis at either observation time point. In contrast, administration of 0.2 mL, even in diluted form, was associated with an increased incidence of vascular occlusion. Overall, these findings indicate that Neofilera(R) presents a lower embolic risk when injected at a volume of 0.1 mL, whereas higher injection volumes may substantially increase the likelihood of embolic complications.

BackgroundDural defects, either from trauma, tumor resection, surgical approaches, genetics, or spontaneously represent a significant clinical challenge in neurosurgery. Given the established efficacy of riboflavin-sensitized ultraviolet-A (UVA) photo-crosslinking in ophthalmology, this study investigated its feasibility and dose-response characteristics as a novel strategy to biochemically augment the mechanical integrity and strength of ex-vivo bovine dura mater. MethodsForty bovine dura mater specimens were treated ex vivo using riboflavin concentrations of 2, 4, or 8 mM combined with UVA irradiation at 0.3 or 3 mW/cm{superscript 2}. PBS-treated specimens exposed to UVA served as controls. Atomic force microscopy nanoindentation was used to measure the local elastic modulus in matched regions before and after treatment, enabling paired assessment of treatment-induced mechanical changes while minimizing inter-sample variability. Post-treatment stiffness, fold-change from baseline, and riboflavin dose-response relationships were analyzed statistically. ResultsBaseline elastic moduli were equivalent across all groups (mean approximately 52 kPa, p=0.92). While UVA alone caused a modest approximately 2- to 3-fold stiffness increase, riboflavin-UVA treatment produced a dramatic, concentration-dependent effect. The highest treatment (8 mM RF, 3 mW/cm{superscript 2} UVA) increased the elastic modulus 150-fold, from approximately 53 kPa to approximately 8,000 kPa. Post-UV stiffness exhibited a strong linear relationship with riboflavin concentration (R{superscript 2} = 0.994), indicating a precisely titratable crosslinking effect. All treatment conditions were statistically distinguishable (p < 0.001). ConclusionRiboflavin-sensitized UVA crosslinking substantially increases the nanomechanical strength of ex vivo bovine dura mater in a controllable, dose-dependent manner. These findings establish a proof of concept for biochemical reinforcement of dural tissue that might be used clinically. As a next step evaluation using human dura, macroscopic biomechanical testing, penetration-depth analysis, and safety assessment is warranted.